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Abstract 



C*~> . We apply an analytic description to the inclusive decay of the r lepton. We 

argue that this method gives not only a self-consistent description of the pro- 
| cess both in the timelike region by using the initial expression for R T and in 

' the spacelike domain by using the analytic properties of the hadronic corre- 

■ lator, but also leads to the fact that theoretical uncertainties associated with 

unknown higher-loop contributions and renormalization scheme dependence 
q ! can be reduced dramatically. 

O 

CL(; I. INTRODUCTION 

The ratio of hadronic to leptonic widths for the inclusive decay of the r-lepton, R T = 
r(r~ — > hadrons u T )/T(r~ — > £ v T \ gives important information about the QCD running 
coupling at relatively small energy scales. The theoretical analysis of the hadronic decay 
of a heavy lepton was performed in [|lj before the experimental discovery of the r-lepton 
in 1975. Since then, the properties of the r have been studied very intensively. Numerous 
publications are devoted to the QCD description of the inclusive decay of the r-lepton and 
determination of the QCD running coupling a s at the r mass scale. A detailed consideration 
of this subject has been given in ||. Recently, an updated QCD analysis has been performed 
by the ALEPH ]3[] and OPAL M collaborations, where applications of different theoretical 
approaches to the r-decay have been analyzed. 

At present, the _R T -ratio is known experimentally to high accuracy, ~ 0.5%. Never- 
theless, the value of a s extracted from the data has a rather large error, in which the- 
oretical uncertainties are dominant. For example, the ALEPH Collaboration result is 
a s (M T = 1.777 GeV) = 0.334 ± 0.007 expt ± 0.021 theor §. It should be emphasized that 
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nonperturbative terms, the values of which are not well known, do not dominate these un- 
certainties, because their contribution is rather small @-§J. The main difficulty is associated 
with the perturbative description. 

The original theoretical expression for the width r(r~ — > hadrons v T ) involves integra- 
tion over small values of timelike momentum [|TJ. The perturbative description with the 
standard running coupling, which has unphysical singularities, becomes ill-defined in this 
region and some additional ansatz has to be applied to get a finite result for the hadronic 
width. To this end, one usually transforms to a contour representation for R T ||, which 
allows one to give meaning to the initial expression and, in principle, perform calculations in 
the framework of perturbative QCD. Assuming the validity of this transformation it is pos- 
sible to present results in the form of a truncated power series with a s (M T ) as the expansion 
parameter There are also other approaches to evaluating the contour integral. The 

Le Diberder and Pich prescription |7j allows one to improve the convergence properties of the 
approximate series and reduce the renormalization scheme (RS) dependence of theoretical 
predictions. The possibility of using different approaches in the perturbative description of 
r-decay leads to an uncertainty in the value of a s (M T ) extracted from the same experimental 
data. Moreover, any perturbative description is based on this contour representation, i.e., 
on the possibility of converting the initial expression involving integration over timelike mo- 
menta into a contour integral in the complex momentum plane. To carry out this transition 
by using Cauchy's theorem requires certain analytic properties of the hadronic correlator 
or of the corresponding Adler function. However, the required analytic properties are not 
automatically maintained in perturbative QCD resummed by the renormalization group. 
It is well known that at the one-loop level the so-called ghost pole occurs in the invariant 
charge. Higher-loop corrections do not solve this problem, but merely add some unphysical 
branch points. The occurrence of incorrect analytic properties in the conventional pertur- 
bative approximation makes it impossible to exploit Cauchy's theorem in this manner and 
therefore prevents rewriting the initial expression for i? T in the form of a contour integral in 
the complex momentum-plane. 

In this paper we will use the analytic approach proposed in || (see also |J for de- 
tails). Being inspired by Kallen-Lehmann analyticity, which is based on general principles 
of quantum field theory, this method ensures that the running coupling possesses the correct 
analytic properties, leads to a self-consistent definition of the effective charge in the time- 
like region [|T0| , |TT| (which cannot be a symmetrical reflection of the spacelike one [Q), and 
provides equality between the initial _R T -expression and the corresponding contour represen- 
tation |TJ]. A distinguishing feature of the analytic approach is the existence of a universal 
infrared limiting value of the analytic running coupling at q 2 = which is independent of 
both the QCD scale parameter A and the choice of renormalization scheme. This limiting 
value is defined by the general structure of the Lagrangian and turns out to be stable with 
respect to higher-loop corrections in contrast to the corresponding quantity in conventional 
perturbation theory (PT). The higher- loop stability of the analytic perturbation theory 
(APT) holds also for physical observables P"5-|iT)| . 

However, it is not sufficient to study the stability with respect to higher-loop corrections; 
one must also investigate the stability with respect to choice of renormalization scheme. This 
is also essential in order to estimate the uncertainty of the results obtained. The theoretical 
ambiguity which is connected with higher-loop corrections and with RS dependence becomes 
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considerable at low energy scales (see, e.g., 
ical version of perturbative QCD 



||17|| ). The APT method, as an invariant analyt- 
improves the situation and gives very stable results 
over a wide range of renormalization schemes. This has been demonstrated for the e + e~ 



annihilation ratio [|14] and for the Bjorken |15| and Gross-Llewellyn Smith [16| deep inelastic 
scattering sum rules. 

The main aim of the paper is a study of the RS dependence which appears in the 
description of the inclusive r decay within the APT approach. We will consider the R T - 
ratio at the next-to-next-to-leading order (NNLO) and the next-to-leading order (NLO) and 
compare results obtained with those of standard perturbation theory. 



II. QCD PARAMETRIZATION OF R T 

The ratio of hadronic to leptonic r-decay widths can be written as 

R T = 3 S EW (\V ud \ 2 + |VU 2 )(l + <*qcd), 



(1) 



where S EW = 1-0194 ± 0.0040 [|IJ is the electroweak factor, \V ud \ = 0.9752 ± 0.0007 and 
\V US \ = 0.2218 ± 0.0016 [20| are the CKM matrix elements, and 5qcd is the QCD correction 
(see for details). 

We first introduce some definitions: Im II ~ 1 + r for the hadronic correlator II (g 2 ) and 
D ~ 1 + d for the Adler function D(q 2 ). Then for massless quarks one can write 5qcd as an 
integral over timelike momentum s: 



JQCD 



M 



t ds ( s V 



1 + 2 W:. 



r(s). 



(2) 



Within the conventional perturbative approximation of r(s) this integral is ill-defined due 
to unphysical singularities of the running coupling lying in the range of integration. The 
most useful trick to rescue the situation is to appeal to analytic properties of the hadronic 
correlator IT(g 2 ). This opens up the possibility of exploiting Cauchy's theorem by rewriting 
the integral in the form of a contour integral in the complex g 2 -plane with the contour being 
a circle of radius M 2 : 



Jqcd 




1 - 



Ml 



1 + m H*> 



(3) 



Starting from the contour representation @ the PT description can be developed in the 
following two ways (see, e.g., plf ). One is Braaten's approach in which the quantity ([|) 
is represented in the form of truncated power series with the expansion parameter a s (M 2 ). 
The NNLO representation for 5qcd is written as follows 



<5qcd = a T + n a 2 + r 2 a 3 T , 



(4) 



where a T = a s (M 2 )/7r. The coefficients r x and r 2 in the MS scheme with three active flavors 
are r\ = 5.2023 and r 2 = 26.366 The running coupling satisfies the renormalization 

group equation: 
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Q (J/CI/ b r\ . rf . 

fJ> -jj^ = I 1 + °i a + °i a )■> ( 5 ) 

where b, c\ and c 2 are the /3-function coefficients. For three active flavors b = 9/2, c\ = 16/9 
and cp = 3863/864. 

In the approach of Le Diberder and Pich (LP) [[/J, the PT expansion is applied to the 
d-function[] 

d(q 2 ) = a(q 2 ) + d ia 2 (q 2 ) + d 2 a 3 (q 2 ) , (6) 
where in the MS-scheme d^ s = 1.6398 and d^ s = 6.3710 for three active quarks. 



Substituting Eq. @ into Eq. (^|) leads to the following expansion, which is not a power 
series in a, 

^qcd = (a) + d x A (2) (a) + d 2 A {3) (a) (7) 

with 

v ; 2ixi J\z\=m? z \ M 2 J \ M 2 J K ' K ' 

As noted above, transition to the contour representation @ requires certain analytic 
properties of the correlator, namely, that it must be an analytic function in the complex 
g 2 -plane with a cut along the positive real axis. The correlator parametrized by the PT run- 
ning coupling does not have this virtue |2"3"| , |T3"|] . Moreover, the conventional renormalization 
group method determines the running coupling in the spacelike region, whereas the initial 
expression (|2|) for R T contains an integration over timelike momentum. Thus, we are in 
need of some method of continuing the running coupling from the spacelike to the timelike 



region that takes into account the proper analytic properties of the running coupling []TT 



Because of this failure of analyticity, Eqs. (0) and (0) are not equivalent in the framework 



of PT [0] and if one remains within PT, nothing can be said about the errors introduced 
by this transition. 

The analytic approach may eliminate these problems. To make our analysis more trans- 
parent and to demonstrate clearly the differences between the consequences of the PT and 
APT methods, we restrict our consideration here to massless NNLO. The NNLO analysis 
can be performed in a more rigorous way without model assumptions that allows us to avoid 
minor details and exhibit the principal features of the APT approach. Thus, other effects, 
such as nonperturbative terms, higher-loop corrections, and renormalon contributions lie 
outside of the purpose of this paper. Note, however, that the NNLO approximation is ade- 
quate to describe the actual physical situation because numerically the corresponding terms 
give the principal contribution to the i? r -ratio. 



1 We use the definition q 2 < in the Euclidean region. We have made a few changes in notation 
from that given in now a = a s /ir, and consequently d\ and d 2 are what we called d 2 and d% 
previously. 
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The function d(q 2 ), which is analytic in the cut g 2 -plane, can be expressed in terms of 
the effective spectral function p(a), the basic quantity in the APT method, 

d{q 2 ) = - -^P(a). (9) 
n Jo a — q A 

The connection between the QCD corrections to the D- and i?-functions can be written 
down in the form of the dispersion integral 

d{q 2 ) = -q 2 r—^-,r(s), (10) 
Jo [s — q l ) 

which is inverted by the following formula 



1 rs+ie rj 7 

rW = ~/ -d(z). (11) 

LTXl Js-ie z 

Here, the contour lies in the region of analyticity of the .D-function. In terms of p{a) the 
function r(s) defined for timelike momenta can be expressed as follows [ TOfl : 



. , 1 r°° da , . , . 

r(s) = - —p{a). (12) 
n Js a 

Eqs. (|9p and (0) determine the QCD corrections d(q 2 ), which is defined in the Euclidean 
(spacelike) region of momenta, and r(s) defined for the Minkowskian (timelike) argument, 
in terms of the spectral function p(a). For <5qcd, using Eq. (0) or equivalently Eq. (|3|), in 
terms of p(cr), we findQ 

1 r 00 da . , 1 f M r da ( a \ 3 / a 



S -—L T^-nk A 1 - Ml) [ 1+ m})^- < 13 > 

In the APT approach, the spectral function is defined as the imaginary part of the 
perturbative approximation to d pt (q 2 ) on the physical cut: 

p{°) = Qo{<r) + d X Qi{a) + d 2 Q2{<?), (14) 

where 

g n (a) = lm[a n+1 (a + ie)} . (15) 

Substituting Eq. flU) into Eq. (0), we can rewrite 5 an in the form of the APT expansion 

5^ = 5® +d l 5^+d 2 5^. (16) 

Note that the APT representations of the <i-function and the QCD correction 5qcd are not 
in the form of power series. 



2 To distinguish APT and PT cases, we will use subscripts "an" and "pt" . 
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The function Qq{(j) in Eq. flUj) defines the analytic spacelike running coupling 



[ /-oo da 

«an(? 2 ) = - / nQo(°)- (17) 

7T Jo o — q A 

In the one-loop approximation it leads to || 

« an (g 2 ) = a pt (g 2 ) + ^^j- ( 18 ) 

Unlike the one-loop PT running coupling, a pt (q 2 ) = 2/61n (— q 2 /A 2 ) , the analytic running 
coupling (|18|) has no unphysical ghost pole and, therefore, possesses the correct analytic 
properties, arising from Kallen-Lehmann analyticity that reflects the general principles of 
the theory. The nonperturbative (non-logarithmic) term, which appears in the analytic 
running coupling, does not change the ultraviolet limit of the theory and thus the APT and 
the PT approaches coincide with each other in the asymptotic region of high energies. 

Thus, the APT approach provides a self-consistent description of the hadronic r decay. 
This description can be equivalently phrased either on the basis of the original expression 
(@), which involves the Minkowskian quantity r(s), or on the contour representation @, 
which involves the Euclidean quantity d(q 2 ). 

An important feature of the APT approach is the fact that d an (q 2 ) and a an (q 2 ) have a 
universal limit at the point q 2 = 0. This limiting value, generally, is independent of both 
the scale parameter A and the order of the loop expansion being considered. Because <i an (0) 
and a an (0) are equal to the reciprocal of the first coefficient of the QCD /5-function, they 
are also RS invariant (we consider only gauge- and mass-independent RSs). The existence 
of this fixed point plays a decisive role in the improved convergence properties relative to 
PT and in the very weak RS dependence of our results. 

To find the analytic function d(q 2 ) involved in Eq. (|3|), we solve the transcendental 
equation for the running coupling 

Hi)-"^ d?S - dl+1 a +Cl,n & +Fm{ah <19) 
where at the NLO 

^m = Wtt^)' (20 » 

and at the NNLO 

F« M = ^M + c 2 /; {l + Cix){l dX +CiX + C2xr (2D 

on the physical cut lying along the positive real axis in the complex g 2 -plane and then use 
Eqs. (|14D, (|15|) and @. Eq. (|19D holds in any MS-like renormalization scheme and allows 
us to normalize the results obtained by using the scale parameter Aj^g. Having found A^g, 
we can study how <5 an varies with a change of renormalization scheme. To do that one has 
to select parameters which determine the RS. The function d(q 2 ) in Eq. (^) is parametrized 
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by a set of RS- dependent parameters. There are RS invariant combinations which constrain 
these parameters |27| . At the NNLO there are two RS-invariant quantities; the first of them 
expresses an energy dependence, the second is just a number 

lo 2 = c 2 + d 2 - c x d x - dx , (22) 

which in our case equals 5.2378. Here, c\ is RS invariant and we can choose d\ and c 2 as 
independent variables, which define some RS. 

There are no fundamental principles upon which one can choose one or another preferable 
RS. Nevertheless, a natural way of studying the RS dependence is to supplement results in a 
certain scheme with an estimate of the variability of the predictions over a range of a priori 
acceptable schemes specified by some criterion. In [^8| it was proposed to consider the class 
of 'natural' RSs, which obey the condition 

\c 2 \ + \d 2 \+c 1 \d l \ + d l 2 <C\uj 2 \. (23) 

This inequality is called the "cancellation index criterion" which means that the degree of 
cancellation in the second RS invariant ( P2"D should not be too large. To define a boundary 
of 'acceptable' schemes which is defined by the value of the cancellation index C, we will 
require no more cancellation than that which occurs in the scheme obeying the principle of 
minimal sensitivity (PMS) [p9[] , which leads to C — 2. 



III. APT: CONVERGENCE PROPERTIES AND RS DEPENDENCE 

For various physical quantities, the APT approach allows one to construct a series that 
has improved convergence properties as compared to a perturbative expansion. To demon- 
strate this fact for the hadronic r-decay, we compare the convergence properties of the 
PT expansions (f|) and (0) on the one hand, and the APT approach given by Eq. (|16"|) 
on the other hand. For our calculation we take as input the TAU'98 conference value: 



R T = 3.642 ±0.019 |2§, which is consistent with the PDG'98 fit R T = 3.642 ±0.024 |20|. In 
Table | we present NNLO results obtained by the methods mentioned above for the central 
experimental value in the MS scheme. The relative contributions of higher-order terms de- 
pends on the method which is applied. The convergence properties of the APT expansion 
seem to be much improved compared to those of the PT expansions. 

The values of the scale parameter A^ and the coupling a s (M 2 ) obtained from above PT 
and APT expansions are noticeably different from each other. The corresponding numerical 
estimations are given in Table PL in which, in order to clarify the situation concerning 
higher-loop stability of different expansions, we also present the NLO result. This table 
demonstrates that the theoretical ambiguity, which associated with different versions of the 
perturbative description, leads to a rather large uncertainty, apx L (Br) — a PT L (LP) = 0.012. 
At the same time the experimental error is Aa expt = 0.007-0.009 The distinction 

between NLO and NLLO running coupling values is 12% for PT (Br) and 5% for PT (LP) 
approaches, while for the APT approach it is less than 0.5%. 

The non-logarithmic terms, which ensure the correct analytic properties and allow a self- 
consistent description of r decay, turn out to be very important for the numerical analysis 
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and influence in an essential way the value of A parameter extracted from the data. Indeed, 
at the one-loop level one can write a simple relation: £ an (A) ~ 5p t p (A) — (2/5)A 2 /M 2 . The 
second term, which is 'invisible' in the perturbative expansion, turns out to be numerically 
important []26| (see the detailed discussion in [|ll|]). Note also that there is a difference 



between the shapes of the analytic and perturbative running couplings, for example, a an (A = 
907 Me V) = 0.403, while at the same scale, the value of the perturbative coupling much 
larger, a pt (A = 907 Me V) = 0.796. Here the question may arise, how is the large APT 
value of A consistent with high energy experimental data? We have estimated the ratio of 
hadronic to leptonic Z-decay widths, Rz, using the above value of A an and the matching 
procedure proposed in [11]. We obtained the value Rz = 20.82, which lies within the range 
of experimental errors; for example, the PDG'98 average is Rz = 20.77±0.07 [pop . This fact 



can be understood if one takes into account that there are differences between the shapes of 
the analytic and perturbative running couplings and also in the terms of the corresponding 
series. 

We found that the value of <5 an depends so slightly on A^g that a 0.9% error in R T gives 
18% error in the value of Aj^g. (This is the reason why the errors in the values of A^jg and 
a(M T 2 ) given by APT are larger than those in PT.) We illustrate this feature in Table FTI . 



According to the table, when we change M 2 /A 2 from 2.0 to 6.5 (corresponding to a variation 
of A from 1.256 GeV to 0.697 GeV), 5 an is only altered by about 20%. The sensitivity to 
A^jg increases as M 2 /A 2 gets smaller. 

Consider now the RS dependence of the APT result and compare it with the perturba- 
tive LP approachFj which of the two PT schemes is more preferable from the viewpoint of 
sensitivity to RS dependence. In the framework of the PT, the RS dependence of <5qcd has 
been discussed in detail in [|!| . 

In the MS scheme we adopt <5 an , P t = 0.1906 and consider some RS belonging to the 
domain described above [see Eq. (p3|) 1- Take two schemes, A and B, located at the two lower 
corners of the boundary of the domain (see Fig. p, i.e., they have the same cancellation 
index as does the PMS scheme, with A = (—1.6183,0) and B = (0.9575,0), where the first 
coordinate is d\ and the second is C2- Then for the PT case in NNLO we get 5 p t{A) = 0.2025 
and Spt(B) = 0.1911. Therefore, even for this sufficiently narrow class of RS the perturbative 
approach gave a 6% deviation in 5qcd that corresponds to a RS uncertainty for the running 
coupling value in the MS-scheme of Aap\ s = 0.0153. The difference between APT results 
is much smaller: 5 an (A) = 0.1890 and <5 an (£>) = 0.1905, and we have only 0.8% deviation, 
which corresponds in the MS-scheme to a RS uncertainty of Act^f = 0.0035. The similar RS 
stability holds also at the two-loop level: one has a 5% deviation in the PT case and only a 
0.4% for the APT. We display our NNLO results in the form of a contour plot, in Fig. |I| 

It is worthwhile to analyze some schemes lying outside the domain considered above 
with the relatively small value of the cancellation index C ~ 2. Among them there is, for 
instance, the commonly used MS scheme which does not belong this domain. In |29|] it was 



shown that the so-called V scheme |30| lies very far from the domain described above and 
gives so a large value of S pt that it cannot be used at this low energy. For the V scheme we 



3 The LP approach is often called the contour-improved fixed-order PT (CIPT) 
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have d\ = —0.109 and C2 = 26.200. The three-loop perturbative result is S pt (V) = 0.3060 
that corresponds to about a 61% deviation from the MS scheme. On the other hand, if we 
turn to APT we have 5 an = 0.1902, i.e., only about a 0.2% deviation from the MS scheme. 
So the V scheme is still useful at this energy in APT. 

In PT at high energies the weak RS dependence is a consequence of the small value of the 
coupling constant. At lower energies the uncertainty increases. In APT, at high energies, 
the situation is the same. However, at low energies the theory has a universal RS-invariant 
infrared limiting value <i an (0), which restricts the RS ambiguity over a very wide range of 
momentum. Another way to illustrate the remarkable stability of APT is to calculate the 
spectral functions g n ((j) given by Eq. (|T5|); one sees that Q\{cr) is much smaller than Qq{o) over 
the whole spectral region. The same statement is true for the relationship between £>i(cr) 
and £2(0") • This monotonically decreasing behavior reduces the RS dependence strongly, 
since the perturbative coefficients d\ and d-i in expression fllil) for p(a) are multiplied by 
these functions. For the MS scheme, this situation is demonstrated in Fig. |2|. 



IV. CONCLUSION 

We have considered inclusive r-decay in three-loop order within analytic perturbation 
theory concentrating on the analysis of theoretical uncertainties coming from the perturba- 
tive short distance part of the QCD correction to the i? T -ratio, which defines the principal 
contribution to this physical quantity. For the low energy r-mass scale, the main source 
of theoretical uncertainties results from the inevitable truncation of the perturbative series, 
which leads to the essential RS dependence and higher loop sensitivity of the theoretical 
predictions. In order to resolve this problem within the conventional perturbative approach 
it is possible to try, in principle, to compute higher loop contributions. However, even if 
this were to be done, one has to keep in mind that from the rigorous point of view it will 
hardly be sufficient because the series is asymptotic, and, in any finite order, the analytic 
properties of the hadronic correlator, which arise from general principles of the theory, are 
violated. Thus, to resolve this problem one has to use a modification of the perturbative 
expansion at low energy scales. 

Here, we have applied the analytic approach which is not inconsistent with the general 
principles of quantum field theory and which opens up the possibility of reducing the the- 
oretical uncertainties associated with short distance contributions mentioned above. Let 
us summarize the important features of this method: (i) the method maintains the cor- 
rect analytic properties and leads to a self-consistent definition of the procedure of analytic 
continuation from the spacelike to the timelike region; (ii) the APT approach has much 
improved convergence properties and turns out to be stable with respect to higher-loop cor- 
rections; (iii) the RS dependence of the results obtained is reduced drastically. For example, 
the V scheme, which gives a very large discrepancy in standard perturbation theory, can 
be used in analytic perturbation theory without any difficulty and the APT predictions are 
practically RS independent over a wide region of RS parameters. 

The nonperturbative power corrections coming from the operator product expansion (in 
this connection see a discussion in J3~T| , |3"2"| | ) , renormalon and other effects are beyond the 



scope of our present consideration. Note, however, that the process of enforcing analyticity 
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modifies the perturbative contributions by incorporating some nonperturbative terms. The 
form of the APT running coupling and the non-power structure of the APT expansion are 
essentially different from the PT ones. Numerically, this difference becomes very important 
in the region less than a few GeV and in order to get the same physical quantity the 
contribution of power corrections should also be changed. 

The value of Aapt is very sensitive to the experimental value of R T . For example, as has 



been demonstrated in [I3| the use of the value of R T obtained by the CLEO collaboration |SB 
gives a value of the scale parameter some 30% smaller than that found here. Note also that 
the renormalon contribution influences the value of A extracted from the r data (see [3?| for a 
review). Within the usual approach, renormalons reduce the value of a s (M^) by about 15%. 
A similar situation holds also in APT and for the nonperturbative a-expansion approach |35| . 
which allows one, as in APT, to maintain the required analyticity [23|| . These two analytic 
approaches often lead to rather similar consequences. For example, they allow one to get 
a good description of experimental data corresponding to the Euclidean and Minkowskian 
characteristics of the process of e + e~ annihilation into hadrons down to the lowest energy 



scale PPl. 



Pure massless APT analysis, which has been performed here, leads to an unusually large 
value of the QCD scale parameter A as compared to the conventional PT value. This 
is connected with the presence of nonperturbative contributions that appear in the APT 
method which have a negative relative sign. The effects mentioned above can change the 
value of the scale parameter extracted from the r data. However, this fact is not relevant for 
the essential conclusion which we have claimed in this paper, that the APT method provides 
predictions which are stable with respect to the choice of renormalization scheme and to the 
inclusion of higher loop corrections. Thus, the analytic approach discussed here is not in 
conflict with the general principles of the theory and allows one to reduce the uncertainties 
of theoretical predictions drastically. 
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FIGURES 




FIG. 1. Contour plot of values of £ an at the three- loop order as a function of RS parameters d\ 
and C2- The dashed line indicates the boundary of the domain, defined by Eq. (|^) with C = 2, 
the heavy points are the positions of the A, B and MS schemes. 
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TABLES 



TABLE I. Successive loop contributions to the PT and APT expansions for 

Rr/[SS EW (\V ud \ 2 + \V us \ 2 )]. 



TVTpthod of dpscrintiori 


Expansion terms 


PT (Br) @ 
PT (LP) § 

A TD r T (loll 

API |§ 


1 + 6% = 1 + 0.104 + 
1 + <^ t p = 1 + 0.148 + 
1 + 5 an = 1 + 0.167 + 


0.056 + 0.030 
0.030 + 0.012 
0.021 + 0.002 


TABLE II. QCD parameters extracted from R T 


= 3.642 ± 0.019 [25] in the MS scheme. 


Approximation Method 


Ams (MeV) 


a(M T 2 ) 


NNLO PT (Br) 

PT (LP) 
APT 


366 ± 14 
391 ± 16 
907 ± 94 


0.328 ± 0.007 
0.340 ± 0.008 
0.403 ±0.015 


NLO PT (Br) 

PT (LP) 
APT 


492 ± 17 
465 ± 19 
954 ± 90 


0.371 ± 0.009 
0.358 ± 0.009 
0.404 ± 0.014 


TABLE III. NLO and NNLO predictions for 5 an in the MS scheme. 


M r 2 /A 2 c e NLO 


m t 2 /a 2 c° 


^NNLO 



2.0 0.2090 0.2106 4.5 0.1820 0.1857 

2.5 0.2016 0.2039 5.0 0.1785 0.1824 

3.0 0.1955 0.1983 5.5 0.1753 0.1795 

3.5 0.1904 0.1935 6.0 0.1724 0.1767 

4.0 0.1859 0.1894 6.5 0.1698 0.1743 
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